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Abstract 

> 

We review the present status of the QCD corrected cross sections and kinematic distributions 
for the production of a Higgs boson in association with bottom quarks at the Fermilab Tevatron 



and CERN Large Hadron Collider. Results are presented for the Minimal Supersymmetric Stan- 
dard Model where, for large tan (3, these production modes can be greatly enhanced compared to 



J3, 

Qh! the Standard Model case. The next-to-leading order QCD results are much less sensitive to the 

<j_) | renormalization and factorization scales than the lowest order results, but have a significant depen- 



dence on the choice of the renormalization scheme for the bottom quark Yukawa coupling. We also 



investigate the uncertainties coming from the Parton Distribution Functions and find that these 

h : 

uncertainties can be comparable to the uncertainties from the remaining scale dependence of the 
next-to-leading order results. We present results separately for the different final states depending 
on the number of bottom quarks identified. 
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I. INTRODUCTION 



One of the most exciting prospects of particle physics is the discovery of the source of 
Electroweak Symmetry Breaking (EWSB). In the Standard Model (SM), a single SU(2)l 
complex scalar doublet is responsible for both gauge boson masses (via the Higgs mechanism) 
and fermion masses (via Yukawa interactions). The only scalar physical remnant of the 
process is the Higgs boson (h) whose mass is a free parameter of the model. In contrast, 
the Minimal Supersymmetric Standard Model (MSSM) contains two complex SU(2)l scalar 
doublets, as required to give masses to both up- and down-type quarks and to cancel gauge 
anomalies. Once EWSB occurs in the MSSM, the scalar sector contains five physical Higgs 
bosons: a light scalar h°, a heavy scalar H°, a pseudoscalar A and a pair of charged bosons 
H ± . Finding experimental evidence for the existence of one or more Higgs bosons is, thus, a 
major goal of current and future collider experiments. In addition, measuring the couplings 
of the Higgs boson(s) to gauge bosons, leptons and quarks would be imperative to uncover 
the true structure of the Higgs sector. 

Direct searches for Higgs bosons at LEP2 have placed lower mass limits on both the 
SM Higgs (M h > 114.4 GeV at 95% c.l.) [jj and the lightest scalar and pseudoscalar Higgs 
bosons of the MSSM (M^o^o^ 93 GeV at 95% c.l.) Q|. In addition, the mass of the lightest 
scalar Higgs boson of the MSSM is also bounded from above by theory to less than ~ 
130 GeV. Finally, precision electroweak measurements, which are sensitive to the effects of 
Higgs boson loop contributions, imply that the mass of the SM Higgs boson should be less 
than ~ 186 GeV at 95% c.l. |3j (~ 219 GeV when including the LEP-2 direct search limit). 
Hence, in both the SM and the MSSM, a Higgs boson should lie in a mass range which can 
be explored at the Large Hadron Collider (LHC) Q]. 

The dominant production mechanism for a SM Higgs boson in hadronic interactions (pp 
or pp) is gluon fusion (gg — > h), which proceeds predominantly through a loop of top quarks. 
Among the production modes with smaller rates, the associated production with electroweak 
gauge bosons (qq — > Vh, where V = W ± , Z°) or with top quark pairs (qq, gg — > tth), as well 
as weak boson fusion (qq — > qqh), play crucial roles. The next-to- leading order (NLO) QCD 
corrections have been calculated for all of the above processes 12. 
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while for gluon fusion and associated Vh production corrections are also known 
at next-to-next-to-leading order (NNLO) in perturbative QCD [17. 18. 19. 20. 21. 22]. 
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FIG. 1: Sample Feynman diagrams for gg — ► bbh and qq — ► production at tree level. 



In the MSSM, the hierarchy of production modes can be significantly different from that 
of the SM. In fact, for large values of tan/? (the ratio of the two MSSM Higgs vacuum 
expectation values), the bottom quark Yukawa couplings to both scalar and pseudoscalar 
Higgs bosons become strongly enhanced and the production of Higgs bosons with bottom 
quarks dominates over all other production modes. 

Recently, the production of Higgs bosons in association with bottom quarks has received 
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24, 251. At tree 



much interest from both the theoretical and experimental communities |2i 
level, the cross section is almost entirely dominated by the sub-process gg — > bbh, with 
only a small contribution from qq — > bbh, at both the Tevatron and LHC (see Fig. |TJ). 
The theoretical prediction of bbh production at hadron colliders, however, involves several 
subtle issues, as will be described in detail in this review, and depends on the number of 
bottom quarks identified, or tagged, in the final state. In case of no or only one tagged 
bottom quark there are two approaches available for calculating the cross sections to bbh 
production, dubbed the four flavor (4FNS) and five flavor (5FNS) number schemes. The 
main difference between these two approaches is that the 4FNS is a fixed-order calculation 
of QCD corrections to the gg, gg-induced bbh production processes, while in the 5FNS the 
leading processes are induced by bg (bg) and bb initial states and large collinear logarithms 
present at tree level as well as at each order of QCD corrections are resummed by using a 
pertubatively defined bottom quark Parton Distribution Function (PDF). The two schemes 
are described in more detail in Section IHBI and III A| respectively. 

For purposes of differentiating the final states, we introduce the following terminology. 
Exclusive bbh production refers to the case when both bottom quark jets are tagged (or are 
at high transverse momentum, p?). Note that, in this case, the computation of the cross 
section solely relies on the 4FNS. Semi-inclusive production (denoted by b(b)h) pertains to 
the case when at least one bottom quark is at high p?. Here, the dominant leading-order 



3 




FIG. 2: Tree level Feynman diagram for bg — * bh in the 5FNS. 



FIG. 3: Tree level Feynman diagram for bb — > h in the 5FNS. 

(LO) process for the 4FNS (5FNS) is gg — ► 6(6) /i (6<? — > 6/i or 6<? — > 6/i, see Fig. |2J). Finally, 
the inclusive production process (denoted by (66) h) involves the case where no bottom 
quark jets are tagged. The dominant LO process for this mode is gg — > (bb)h (66 — > h, 
see Fig. |3J) in the 4FNS (5FNS). It should be mentioned that, although the cross section 
for the inclusive production mode dominates over the less inclusive modes by one to two 
orders of magnitude, tagging bottom quark jets greatly reduces backgrounds, making the 
less inclusive modes experimentally more appealing. Requiring one or two high-pr bottom 
quark jets in the final state also ensures that the detected Higgs boson has been radiated 
off a bottom or anti-bottom quark and the corresponding cross section is unambiguously 
proportional to the bottom quark Yukawa coupling. 

QCD corrected cross sections are available for all three final states and the inclusive and 
semi-inclusive production processes have been computed in both the 4FNS and 5FNS. For 



the inclusive case, the NLO QCD corrected 4FNS |23|, |2(j and the NNLO QCD corrected 
5FNS cross sections have been compared and are found to be in good agreement within 



theoretica 
4FNS 



uncertainties. The NLO predictions of the semi-inclusive cross sections for the 
^ and 5FNS Q nave also been ex t en S ive. y _d and the a„t 



2f 

between the two is spectacular. The compatibility of these two seemingly different calcula- 
tional schemes in the prediction of inclusive and semi-inclusive Higgs boson production rates 
is indeed a powerful check of the theory. The 4FNS and 5FNS represent different perturba- 
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tive expansions of the same physical observables, and therefore should agree at sufficiently 
high order. Indeed, we see that considering the first (second) order of corrections already 
brings agreement between the two schemes within the respective theoretical uncertainties. 
Finally, two independent calculations of the NLO QCD corrections for the exclusive mode 



23, 



26 
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have been compared and agreement has been found 

The above discussion for the production of a scalar Higgs boson with bottom quarks 
applies equally well to the production of a pseudoscalar Higgs boson. In fact, for massless 
bottom quarks, the predictions for bbA° would be identical to those of bbh°/H° given a 
rescaling of the Yukawa couplings. For massive bottom quarks, however, the situation 
becomes complicated due to the 75 appearing in the bbA° Yukawa coupling. The 75 Dirac 
matrix is intrinsically a four- dimensional object and care must be taken in its treatment 
when dimensionally regularizing the calculation. In any case, bottom quark mass effects are 

2 _ _ 

expected to be small, i.e. O(jti), and predictions for bbh°/H° are good indicators for bbA° 

h 

production (modulus Yukawa couplings). 

The remainder of this paper is organized as follows. In Section |Hj we discuss the frame- 
work of the calculations in the four- and five-flavor-number schemes in detail. Given the 
importance it plays in the prediction of the cross section through the Yukawa coupling, we 
also discuss different renormalization schemes of the bottom quark mass. Section ITTT1 details 

□ in 

the SUSY-corrected Yukawa couplings obtained from the program FeynHiggs 30, 31, 32 
that were used in this study to produce results for MSSM Higgs bosons. In Section HV| we 
review the results for inclusive, semi-inclusive and exclusive Higgs boson production with 
bottom quarks. In Section we investigate the uncertainties arising from the PDFs for 
the semi-inclusive bbh production process. These uncertainties are calculated using the al- 



gorithm developed by the CTEQ collaboration 
Section Ell 



33j. Finally, we summarize and conclude in 



II. THEORETICAL FRAMEWORK 
A. Five Flavor Number Scheme 

Final state bottom quarks are identified imposing cuts on their transverse momentum 
(pT>PT in ) and pseudorapidity (\rj\ <ri max ). On the other hand, when a final state bottom 
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quark is not identified, as in the inclusive and semi-inclusive bbh production processes, the 
corresponding integration over its phase space may give rise to large logarithms of the form: 



where m& and \x represent the lower and upper bounds on the integration over the trans- 
verse momentum, p?, of the final state bottom quark. This happens when the final state 
bottom quark is directly originating from the g — > bb splitting of an initial state gluon, 
and corresponds to the collinear singularity that would be present in the case of a gluon 
splitting into massless quarks. The scale \i is typically of 0(M h ) and therefore, due to the 
smallness of the bottom quark mass, these logarithms can be quite large. Additionally, the 
same logarithms appear at every order in the perturbative expansion of the cross section in 
the strong coupling, a s , due to recursive gluon emission from internal bottom quark lines. 
These logarithms could severely hinder the convergence of the perturbative expansion of the 
total and differential cross sections. In the 5FNS the convergence is improved by introducing 
a perturbatively defined bottom quark PDF ^J, [3 S| , defined at lowest order in a s as 



where g(y,fi) is the (experimentally-measured) gluon PDF and P qg is the Altarelli-Parisi 
splitting function for g — > qq given by 

P qg {z) = l -^ + {l-zf] . (3) 



Subsequently, the A b logarithms are resummed through the DGLAP evolution equation 



38 



37 



39|, such that contributions proportional to (a s Ai,) n can be absorbed, to all orders in n, 



into a leading logarithms bottom quark PDF, while subleading logarithms are recursively 
resummed when higher order corrections are considered in the DGLAP equation. Since 
a s Ab is not a small expansion parameter, the use of a bottom quark PDF should improve 
the stability of the total and differential cross sections for inclusive and semi-inclusive Higgs 
boson production with bottom quarks. 

The 5FNS is based on the approximation that any spectator bottom quarks (i.e., bottom 
quarks which are not tagged in the final state) are at small transverse momentum, since this 
is the region where the A& logarithms dominate over other ^-dependent contributions. At 
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lowest order, the spectator quarks are produced with zero pr, and a transverse momentum 



34, 



35]. 



spectrum for the outgoing bottom quarks is generated at higher orders 

With the use of a bottom quark PDF, the 5FNS effectively reorders the perturbative 
expansion to be one in a s and A^ 1 . To see how this works, let us consider the perturbative 
expansion of the inclusive process bb — > h (Fig. |SJ), which is at tree level of order a^A|. 
At NLO, the virtual and real corrections to the tree level process make contributions of 
0(a^Al). However, at NLO, we must also consider the contribution from bg — > bh where 
the final state bottom quark is at high p?. This process makes a contribution of order 
ctgAfc and is, thus, a correction of 0{K b 1 ) to the tree level cross section. Similarly, at NNLO, 
among the myriad of radiative corrections of 0{a i s A 2 ), we must also include the contribution 
from the process gg — > bbh, where both b and b are at high p?. The contribution from these 
diagrams are of order a 2 s , and are, thus, 0(A^ 2 ) (or NNLO) corrections to the tree level 



process bb^h [36, 40]. 

The above discussion for bb — > h also applies to the perturbative expansion of bg — ► bh. 
In this case, the tree level process is of order a? A}, and the contribution from gg — > bbh is a 
NLO correction of 0(A^) y. 

Finally, we observe that, in the existing 5FNS calculations, all bottom quark masses other 
than the one appearing in the Yukawa coupling are set to zero. 



B. Four Flavor Number Scheme 



In the 4FNS, the initial state quarks are constrained to be the four lightest quarks, i.e. 
there are no bottom quarks in the initial state. No kinematic approximations are made and 
the cross section for pp(pp) — > bbh is computed at fixed order in QCD without resumming 
higher-order collinear logarithms. Moreover, the bottom quark mass is always considered to 
be non-zero. Similarly to the 5FNS, for the less inclusive final states discussed above (b(b)h 
or bbh), identification cuts are imposed on the final state bottom quark(s) which constrain 
its transverse momentum and pseudorapidity. Applying these cuts to a final state bottom 
quark eliminates large logarithms that may appear from gluon splitting and it ensures that 
the bottom quark jets can be tagged experimentally. While the inclusive and semi-inclusive 
production of a Higgs boson with bottom quarks can be calculated in both the 4FNS and 
5FNS, the exclusive production can only be calculated in the 4FNS, since both final state 
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bottom quarks are tagged. 

The NLO QCD corrections to the hadronic process pp(pp) — > bbh in the 4FNS consist of 
calculating the Q (ou ) virtual and real QCD corrections to the tree level (partonic) processes 
29[. We note that the NLO calculation of p p( pp) — » bbh is identical to the 



QQ,99 



bbh 



26 



calculation of the NLO QCD corrections to pp(pp) — > tth jll, 
interchange of the top quark and bottom quark masses (m t ■ 



12 



13 
m b ). 



with the global 



C. Definition of the bottom quark mass 

One potential source of theoretical uncertainty in the calculation of pp(pp) — ► bbh involves 
the renormalization of the bottom quark mass. The bottom quark mass counterterm has 
to be used twice in the renormalization of the calculation: once to renormalize the bottom 
quark mass appearing in internal propagators and once to renormalize the SM bottom 
quark Yukawa coupling gr^j-f = rrib/v, where v = (v^GV) -1 / 2 = 246 GeV. Indeed, if one 
only considers QCD corrections, the counterterm for the bottom quark Yukawa coupling 
coincides with the counterterm for the bottom quark mass, since the vacuum expectation 
value v is not renormalized at 1-loop in QCD. 

The renormalization scheme dependence has been studied in Ref. j^, in the 4FNS, for 
the case of the exclusive bbh production. Two schemes for the renormalization of the bottom 
quark mass have been chosen: the MS scheme and an on-shell (OS) scheme. The two are 
perturbatively consistent at NLO with the only difference being at higher orders and, thus, 
part of the theoretical uncertainty of the NLO calculation. The main reason for investigating 
this renormalization scheme dependence is the large sensitivity of the MS bottom quark 
mass and the prominent role it plays in the bbh production cross section through the overall 
bottom quark Yukawa coupling. A summary of the results of this investigation are presented 
in Section WU 

The benefit of using the MS scheme for the Yukawa coupling is that it potentially gives 
control over higher-order corrections beyond the 1-loop corrections. This is often reflected 
in the better perturbative behavior of physical observables calculated using the MS bottom 
quark Yukawa coupling. Therefore, unless otherwise stated all the results presented in this 
review are given using the MS scheme for the bottom quark mass and Yukawa coupling. 
This is accomplished by replacing the MS bottom quark mass in the Yukawa coupling with 



S 



the corresponding 1-loop and 2-loop renormalization group improved MS masses: 
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where we take the bottom quark pole mass to be m b = A.62 GeV and 
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are the one and two loop coefficients of the QCD (3 function and mass anomalous dimension 
7 m , while N = 3 is the number of colors and n// = 5 is the number of light flavors. 



III. RADIATIVELY-CORRECTED MSSM BOTTOM YUKAWA COUPLING 

So far, most of the results for bbh production have been obtained in the SM, in order 
to simplify the comparison of different calculations and approaches (for a review see, e.g., 
Ref. H)- 

As stated earlier, though, these processes are relevant discovery modes at hadron 
colliders only for SUSY models with large tan (3. Therefore, we will present most of the 
following results for the case of the MSSM with large tan f3. The prediction for the MSSM 
case can easily be derived from the SM results by rescaling the Yukawa couplings. In the 
4FNS, due to the existence of diagrams where the Higgs is radiated off a closed top quark 
loop, some care must be taken when performing this rescaling. In this review, all 4FNS 
results for the neutral MSSM Higgs bosons, are obtained from the SM results by carefully 
rescaling the contributions with (a SM ) and without top quark loops (a SM — <r| M ) as: 

(_ MSSM \ 2 / MSSM _ A4SSM \ 

Ubbh \ (- t \ _i | Ugh y bb h \ _t ( 7 \ 

n SM J l°SJVl &SM) "T I SM n SM J SM ■ {') 

ybbh / v ytthy bb h / 

Examples of the closed top quark loop contributions described by cr* are shown in Fig. HI 
Note that in the 5FNS, as implemented in the program MCFM |4l| . the corresponding 
contributions to cr* M are considered to be zero. While this could affect the comparison in 
the SM, in SUSY models with large bottom quark Yukawa couplings (and consequently 
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FIG. 4: Sample Feynman diagrams for the closed top quark loop contribution to gg — ► bbh. 
The circled cross denotes all possible insertion of the final state Higgs boson leg, each insertion 
corresponding to a different diagram. 



small top quark Yukawa couplings) the contributions of the closed top quark diagrams is 
very small, and <7* M does not play a significant role. 

In the MSSM the bottom and top quark Yukawa couplings to the scalar, neutral Higgs 
bosons are given by: 



n MSSM 

9bbh° 



n MSSM 

ybbH° 



sin a 



cos (3 



9bbh 



cos a 



9 bbh 



n MSSM 

ym 



r .MSSM 

9itH° 



cos a 



sin (3 



sin a 



9'tth 



9tth 



(8) 



cos(3 -bbn > w gin/3 . 

where g^' h and g"^ are the SM bottom and top quark Yukawa couplings, h° and H° are 
the lighter and heavier neutral scalars of the MSSM, and a is the angle which diagonalizes 

rn 

the neutral scalar Higgs mass matrix J42j. The dominant SUSY radiative corrections to 
bbh production can be taken into account by including the SUSY corrections to the bbh 
vertex only, i.e. by replacing the tree level Yukawa couplings by the radiative corrected 
ones. We follow the treatment of the program FeynHiggs [43] and take into account the 
leading, tan/3 enhanced, radiative corrections that are generated by gluino-sbottom and 
chargino-stop loops, as follows [44]: 

1 



MSSM 

9bbh° 



MSSM 
9bbH° 



~9bbh 



9bbh 



1 + A 6 
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1 + A h 



sin a 

cos (3 
cos a 
cos (3 



A, 



+ A, 



cos a 
sin [3 
sin a 
'sin (3 
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with 



//tan/3 



2a s (m t ) 
3tt 



M~ g I(m~ bi ,m~ b2 , 



m~ g ) + 
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A t I( mil ,m i2 ,fi) 



(10) 
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tan /3 = 40 


M h0 [GeV] 


100 110 


120 


130 


a [rad] 


-1.5063 -1.4716 


-1.3798 


-0.7150 


MSSM 1 SM 

y bbh° 1 y bbh 


33.913 33.823 


33.387 


22.390 


n MSSM 1 SM 

Urn I9tth 


0.0645 0.0991 


0.1899 


0.7553 



TABLE I: Ratios of the MSSM and SM bottom and top quark Yukawa couplings to a light neutral 
scalar, along with the corresponding values for a as computed by FeynHiggs 4^|. The genuine 
SUSY input parameters are chosen as follows: all soft-SUSY breaking masses for squark doublets 
and singlets are set to M SUSY = 1 TeV, Mg = l TeV, A b = A t = 2 TeV, and /j = M 2 = 200 GeV. For 
this choice we find A& = 0.178. 

where h t = rn t (m t )/ y/vf + v\ ~ m t (m t )/ 17 4 GeV, mi ,mf 12 ,Mp denote the sbottom, stop, 
and gluino masses, respectively, A t is the stop and [i the Higgs mixing parameter. The strong 
coupling constant and the running top quark mass, m t , are evaluated at the on-shell top 
mass, m t = 174 GeV. The vertex function I is given by 

[a 2 6 2 ln(^)+6 2 c 2 ln(^) + a 2 c 2 ln(4)] , > 

T a b c) = - 6 c a 11) 

In Tables HI and ITT1 we provide the ratios of the SM and MSSM bottom and top quark Yukawa 
couplings at tan/3 = 40, as provided by FeynHiggs, where we assumed CP conserving 
couplings 1 . In the following, we will use these ratios to derive the MSSM cross sections 
from the SM results as described in Eq. ©. Using FeynHiggs, the MSSM Higgs boson 
masses and the mixing angle a have been computed up to two-loop order. As expected, for 
this choice of MSSM input parameters the bottom Yukawa coupling is strongly enhanced 
and the top Yukawa coupling is suppressed, resulting in a MSSM bbh cross section that is 
about three orders of magnitude larger than the SM cross section. 

IV. NUMERICAL RESULTS 

If not stated otherwise, the numerical results are obtained using CTEQ6M parton distri- 
bution functions for the calculation of the NLO cross section, and CTEQ6L parton distribu- 



1 MSSM Yukawa couplings are also provided by the programs HDECAY and CPSuperH 0E 

11 



tan (3 = 40 



M H0 [GeV] 120 200 300 400 500 600 800 
a [rad] -0.3315 -0.0454 -0.0318 -0.0282 -0.0266 -0.0258 -0.0251 

dblHo/dbih 25 - 787 27 - 338 27 - 356 27 - 360 27 - 362 27 - 363 27 - 364 
9&p/9tfh -°- 3255 -0-0454 -0.0318 -0.0282 -0.0266 -0.0258 -0.0251 

TABLE II: Ratios of the MSSM and SM bottom and top quark Yukawa couplings to a heavy neutral 
scalar, along with the corresponding values for a as computed by FeynHiggs 4^|. The genuine 
SUSY input parameters are chosen as follows: all soft-SUSY breaking masses for squark doublets 
and singlets are set to M SU sy = 1 TeV, M s = l TeV, A b = A t = 25 GeV, and fi = M 2 = 1000 GeV. 
For this choice we find A& = 0.461. 



tion functions for the calculation of the lowest order cross sections |48[. The NLO (LO) cross 
section is evaluated using the 2-loop (1-loop) evolution of a 8 (fj) with af LO (M z ) = 0.118. 
For the exclusive and semi-inclusive channels (bbh and bh + bh production), it is required 
that the final state bottom quarks have pr > 20 GeV and pseudorapidity \rj\< 2.0 for the 
Tevatron and | rj |<2.5 for the LHC. In the NLO real gluon emission contributions, the final 
state gluon and bottom quarks are considered as separate particles only if their separation 
in the pseudorapidity-azimuthal angle plane, AR = a/ (Arj) 2 + (A</>) 2 , is larger than 0.4. For 
smaller values of AR, the four momentum vectors of the two particles are combined into an 
effective bottom/anti-bottom quark momentum four- vector. 

In the following numerical discussion, the SM NLO QCD corrected cross sections for 
semi-inclusive and exclusive bbh production in the 4FNS are taken from Refs. and Q, 
respectively. The 5FNS SM results for semi- inclusive b(b)h production have been produced 
with the program MCFM j^j and for the inclusive case have been taken from Ref. 



A. Results for inclusive (bb) h production 

If the outgoing bottom quarks cannot be observed then the dominant MSSM Higgs pro- 
duction process at large tan/5 is gg — * (bb)h (4FNS) or bb — > h (5FNS). At the LHC, this 
process can be identified by the decays to and t + t~ for the heavy Higgs bosons, H° 

and A , of the MSSM where the Higgs couplings are enhanced at large tan/?. Recently, also 
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0-nlo [pb] Tevatron, y/s = 1.96 TeV o"nlo [pb] LHC, ^ = 14 TeV 




100 105 110 115 120 125 130 150 200 250 300 350 400 



M h o [GeV] M H o [GeV] 

FIG. 5: Total cross sections for pp,pp — > (bb)h (h = h°,H°) in the MSSM with no bottom quark 

jet identified in the final state in the 4FNS (at NLO) and 5FNS (at NLO and NNLO) as a function 

of the light and heavy MSSM Higgs boson masses, at both the Tevatron and the LHC. The error 

bands have been obtained by varying the renormalization and factorization scales as described in 

the text. 



at the Tevatron this process, with h — > r + r , has been used to search for t 



he neutral MSSM 



3, 



and has only 



Higgs boson |25j]. In the 5FNS this process has been computed to NNLO 
a small renormalization and factorization scale dependence. In the 4FNS, the production 
processes cm, qq — ► (bb)h, where the outgoing bottom quarks are not observed, are known at 
NLO 0,13 an d have been recalculated for the purpose of this review. A comparison of the 
total production rates in the 4FNS and 5FNS within the MSSM for tan /3 = 40 is shown in 
Fig. 03 The bands represent the theoretical uncertainty due to the residual scale dependence 
at NLO (NNLO). In the 4FNS they have been obtained by varying the renormalization (/i r ) 
and factorization (fif) scales independently from /^o/4 to /io, where fio = mb + In the 

5FNS, the renormalization scale is fixed to fi r = and the factorization scale is varied 
around the central scale /i/ = M^/4 from O.IM^ to 0.7M/,. There is good agreement between 
the results of the two schemes within their respective scale uncertainties, although the 5FNS 
at NNLO is slightly higher than the 4FNS for all values of the Higgs boson masses. 

B. Results for semi-inclusive b(b)h production 



If a single bottom quark is tagged then the final state is bh or bh. A recent Tevatron 



study 



24j used the search for neutral MSSM Higgs bosons in events with three bottom 
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FIG. 6: Total LO and NLO cross sections in the SM for pp,pp — > b(b)h production in the 4FNS as 
a function of /i = fi r = for M/j = 120 GeV, at both the Tevatron and the LHC. 



quarks in the final state (bh + bh production with h — > bb) to impose limits on the tan /3 
and Af^o parameter space. The effect of the NLO corrections in the 4FNS are illustrated 
in Fig. |H1 where we show the SM results for the LO and NLO total production rates (27 1. 
We see a strong reduction of the scale dependence at NLO. In Fig. [7|we compare the rates 
obtained in the MSSM for tan/3 = 40 for the single 6-tagged events in the 4FNS and 5FNS 
schemes, where we varied separately the renormalization and factorization scales from 0. 2un 
to /io, with /i = rrib + Mh/2. These are obtained by rescaling the SM results of Ref. |27| 
according to Eq. (J7J) to the MSSM scenario discussed in this review. Note that the resulting 
bands only give an estimate of the theoretical uncertainty due to residual scale dependence, 
but do not consider other uncertainties such as PDF uncertainties. A discussion of the 
PDF uncertainties can be found in Section As demonstrated in Fig. [7| the theoretical 
predictions in the 4FNS and 5FNS are fully compatible and, thus, either calculation can 
be used in the experimental analyses. The theoretical uncertainty in both schemes is still 
dominated by the residual scale dependence and the PDF uncertainty. 

In comparing the four and five flavor number schemes it is particularly interesting to 
compare the kinematic distributions. In Figs. |HK1 we compare the results for the transverse 
momentum and pseudorapidity distributions of the light and heavy MSSM Higgs bosons, h° 
and if , in both the 4FNS and 5FNS, at the Tevatron and the LHC. We see, in general, good 
agreement between the two schemes within their respective theoretical uncertainties, except 
in regions of kinematic boundaries. This is particularly dramatic in the y\ distributions in 
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o-nlo [pb] Tevatron, sfs = 1.96 TeV CT NLO [pb] LHC, y/a = lA TeV 




100 105 110 115 120 125 130 150 200 250 300 350 400 450 500 



M h o [GeV] M H o [GeV] 

FIG. 7: Total NLO cross section in the MSSM for pp,pp — > b(b)h production at the Tevatron and 

the LHC as a function of M^o jjo. We varied \i r and Hf independently to obtain the uncertainty 

bands, as explained in the text. The solid curves correspond to the 4FNS, the dashed curves to 

the 5FNS. 



the 5FNS where, around p\ — 20 GeV, a kinematic threshold induced in bg — > bh by the 
cut on the p? of the bottom quark causes the 5FNS NLO calculation to be unreliable. This 
instability can be reabsorbed by using a larger bin size (see inlays), and could therefore be 
interpreted as a sort of theoretical resolution for the 5FNS. The instabilities could be removed 
by a systematic resummation of threshold corrections 
MCFM. 



wh ch is not implemented in 



C. Results for exclusive bbh production 

Finally, we discuss the fully exclusive LO and NLO cross sections for bbh production, 
where both the outgoing b and b quarks are identified. The results of Ref. 29| that are 
presented here have been obtained by using the CTEQ5 set of PDFs. We have checked 
that the conclusions of the numerical discussion presented here do not change when using 
the CTEQ6 set of PDFs. In Fig. EH we show, for M h = 120 GeV, the dependence of the 
LO and NLO total cross sections, calculated in the SM and in the 4FNS, on the arbitrary 
renormalization/factorization scale \x (with fi r — fif — fi). The curves labelled "OS" ("MS") 



use the on-shell (MS) scheme for the bottom Yukawa coupling. The NLO results have 



significantly less sensitivity to the scale choice, both in the OS and in the MS schemes. 
However, while in the OS scheme the plateau region, or region of least sensitivity to the 
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FIG. 8: do/dp\ in the MSSM at the Tevatron and the LHC for M ft o iH o = 120 GeV and i i r = i x F = 
/xo/2 for single 6-tag events. We show the NLO results in the 4FNS (solid) and 5FNS (dashed), 
using two different bin sizes, 2 GeV (left) and 12 GeV (right). 



renormalization/factorization scale, is around /i = /io, in the MS scheme it is shifted towards 
[i = Ho/2, (/i = rrib + Mh/2). It is interesting to note, that the scale choice \i = Mh/4 is 
supported by theoretical studies |3] and corresponds to the point where the NM 
for inclusive Higgs production from bottom quarks is the same as the NLO rate 
Fig. El shows how the MS scheme gives a perturbative cross section that is better behaved 
over a broader range of scales, and therefore preferable. Conservatively, one could interpret 
the difference between the two plateau regions as an additional 10-20% theoretical error. 
In Fig. ^2 we compare the production of the SM Higgs boson with that of the neutral 
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FIG. 9: da/dr/h in the MSSM at the Tevatron and the LHC for M h o jjo = 120 GeV and [i T = f = /io/2 
for single 6-tag events. We show the NLO results in the 4FNS (solid) and 5FNS (dashed). 
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FIG. 10: Total LO and NLO cross sections in the SM for the two 6-tag process, pp,pp — ► bbh, in the 
4FNS as a function of /x = fi r = [J-f for = 120 GeV, at both the Tevatron and the LHC. Results 



are shown for both the MS and OS renormalization scheme for the bottom Yukawa coupling (see 
discussion in Section Pi C|) . 



scalar Higgs bosons of the MSSM and again observe a significant enhancement of the rate 
in the MSSM for large tan/3. The MSSM results have been obtained from the SM results 
as described in Eq. EJby using the rescaling factors of Tables HI HT1 for tan/5 = 40. 

Finally, in Figs. El El we illustrate the impact of NLO QCD corrections on the trans- 
verse momentum and pseudorapidity distribution of the SM Higgs boson and the bottom 
quark by showing the relative correction, da^io I d&LO ~ 1 (4 n percent). For the renormaliza- 
t ion/factorization scale we choose [i = 2fi at the Tevatron and fi = 4/i at the LHC. These 
two scale choices are well within the plateau regions where the NLO cross sections vary the 
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FIG. 11: Total NLO cross sections in the SM and MSSM for the two 6-tag process, pp,pp —* bbh, 
in the 4FNS for the central scale \i r = [it = fio- Shown is the dependence on the light and heavy 
MSSM Higgs boson masses at the Tevatron and the LHC, respectively. 
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FIG. 12: The relative corrections da NLO /da LO — 1 for the px distribution of the bottom or anti- 
bottom quark with the largest pr (pT,max) (left) and of the SM Higgs boson (pj,) (right) to bbh 
production in the SM at the Tevatron (with y/s = 2 TeV and \i = 2/j.q) and the LHC (with 
s/s = 14 TeV and ^ = 4^ )- 



least with the value of fi. As can be seen, the NLO QCD corrections can considerably affect 
the shape of kinematic distributions, and their effect cannot be obtained from simply rescal- 
ing the LO distributions with a K-factor of cr NLO /cr LO — 1.38 ± 0.02 (Tevatron, /i = 2/i ) and 
Cnlo/vlo = 1-11 ± 0.03 (LHC, /i = 4/i ). The errors only include the statistical uncertainty 
from the Monte Carlo integration. The kinematic distributions have been calculated within 
the SM and using the OS scheme, but we see a similar behavior when using the MS bottom 
quark Yukawa coupling or consider the MSSM case. 
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FIG. 13: The relative corrections da NLO /da LO — 1 for the rj distribution of the bottom quark r#, 
(left) and of the SM Higgs boson (rjh) (right) to bbh production in the SM at the Tevatron (with 
y/s = 2 TeV and u = 2{iq) and the LHC (with -^/i=14 TeV and /z = 4/io)- 



V. PDF UNCERTAINTIES 

Besides the residual renormalization/factorization scale dependence after the NLO correc- 
tions have been included, another major source of theoretical uncertainty for cross section 
predictions at hadron colliders comes from the Parton Distribution Functions. Unfortu- 
nately, PDFs are plagued by uncertainties either from the non-perturbative starting dis- 
tributions used to fit the data or from the perturbative DGLAP evolution to the higher 



energies relevant at hadron colliders |37|, 



38, 



3< 



Recently, several collaborations have introduced new schemes which allow an estimate of 
theoretical uncertainties on physical observables due to the uncertainty in the PDFs. Here, 
we focus on the scheme introduced by the CTEQ collaboration based on the Hessian matrix 
method ||3j] and study the uncertainties of semi-inclusive total bbh production rates that 
are induced by the uncertainties in the PDFs. The details of this method are beyond the 
scope of this review, however, we give a brief explanation below. First, the nominal set of 
PDFs (e.g. CTEQ6) is constructed by fitting a non-perturbative core equation to data from 
low-energy experiments designed to measure PDFs. The core equation, in the method used 
by CTEQ, is parameterized by 20 independent parameters which are dialed to fit the data. 
Once the nominal set is fixed, the 20 parameters are then varied in a well-defined manner 
to produce an additional 40 sets of PDFs. These sets serve as a map of the neighborhood 
around the nominal fit to the data. Indeed, one can then use the 40 sets to estimate the 
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uncertainty from the PDFs on a physical observable in the following way [53] 2 : first, the 
central value cross section Go is calculated using the global minimum PDF (i.e. CTEQ6M). 
The calculation of the cross section is then performed with the additional 40 sets of PDFs to 
produce 40 different predictions, (7j. For each of these, the deviation from the central value is 
calculated to be Aaf = | cr^ — cr 1 when g^Gq. Finally, to obtain the uncertainties due to the 



PDFs the deviations are summed quadratically as Acr ± = y ^ Aof" 2 and the cross section 
including the theoretical uncertainties arising from the PDFs is quoted as o'oHa' 7 - • Recently, 
similar analyses have been performed for the dominant SM Higgs production modes Q and 



inclusive Higgs production modes in both the SM {gg — > h) and the MSSM 
1 . In Fig s. IT41 and ^] we plot the total SM NLO cross section for bg — > bh 

C 



for the ful 
(66 -> h) 

obtained with MCFM [4JJ at the Tevatron and LHC respectively. Here, we compare the 
spread uncertainties from residual scale dependence and the PDFs both for the total cross 
section (top) and the total cross section normalized to the central value calculated with 
CTEQ6M (bottom). 

From Fig. El one can see that, at the LHC, the theoretical uncertainty is dominated by 
the residual scale dependence. Due to the large center of mass (c.o.m.) energy of the LHC, 
the gluons and bottom quarks in the initial state have small momentum fraction (x) values 
and, hence, small PDF uncertainties typically in the 5-10% range. 

In contrast, due to the smaller c.o.m. energy, the PDF uncertainties at the Tevatron 
(Fig. ITlj) are comparable and even larger than the uncertainties due to residual scale depen- 
dence over the full Higgs mass range. The smaller c.o.m. energy results in higher- £ gluons 
and bottom quarks in the initial state which corresponds to large PDF uncertainties in the 
10-30% range. 

Finally, in Fig. we plot the normalized total SM NLO cross sections of bg — > bh and 
gg — > 6(6) h and compare their respective uncertainties due to the PDFs. We see that, at 
both the Tevatron and the LHC, the PDF uncertainties are almost identical for both the gg 
and bg initial states. 



2 We have also performed this analysis using the PDF sets of the MRST collaboration [^j. These sets are 
made up of 30 sets in addition to the nominal fit and, hence, map less of the neighborhood around the 
global minimum. This results in smaller spread uncertainties than the CTEQ analysis. Therefore, we 
only show results using the CTEQ sets and quote these results as an upper limit of the uncertainty from 
PDFs. 
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FIG. 14: Comparison between theoretical uncertainties due to scale dependence and uncertainties 
arising from the PDFs at the Tevatron. In the bottom plot, both uncertainty bands have been 
normalized to the central value of the total cross section ctq. 



VI. CONCLUSIONS 

We discussed the status of theoretical predictions for Higgs boson production in associ- 
ation with bottom quarks for the cases where zero, one or both outgoing bottom quarks 
are identified. These processes are important Higgs boson discovery modes in models with 
enhanced bottom quark Yukawa couplings, such as the MSSM at large values of tan (3. We 
presented results for the MSSM bbh (h = h°, H°) total and differential production rates at 
tan /3 = 40 and discussed the theoretical uncertainties of the state-of-the-art predictions due 
to the residual renormalization/factorization scale dependence at NLO (NNLO) in QCD and 
the uncertainty in the PDFs. We have shown that the total cross section for the inclusive 
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FIG. 15: Comparison between theoretical uncertainties due to scale dependence and uncertain- 
ties arising from the PDFs at the LHC. In the bottom plot, both uncertainty bands have been 
normalized to the central value of the total cross section gq. 

case, pp,pp — > (bb)h, and both the total and differential cross sections for the semi- inclusive 
case, pp,pp — > b(b)h, within a 4FNS and 5FNS are fully compatible within the existing the- 
oretical errors. All important Higgs production processes at hadron colliders are now known 
at NLO, and some even at NNLO, and the theoretical uncertainties of the QCD predictions 
are well under control. 
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